Spin and charge dynamics of the 2D t-J model at intermediate electron densities: 

absence of spin-charge separation 
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We present an exact diagonalization study of the dynamical spin and density correlation func- 
tions in small clusters of t — J model, focussing on the regime of intermediate and low electron 
densities, p e <0.5. In 2D both correlation functions agree remarkably well with the convolution of 
the single-particle spectral function, i.e. the simplest estimate possible within a Fermi liquid pic- 
ture. Deviations from the convolution are shown to originate from symmetry-related selection rules, 
which are unaccounted for in the convolution estimate. For all fillings under consideration, we show 
that the low energy peaks originate from particle-hole excitations between the Fermi momenta, as 
expected for a Fermi liquid. We contrast this with the behaviour in ID, where spin and density 
correlation function show the differences characteristic of spin-charge separation and where neither 
correlation function is approximated well by the convolution. 
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A much discussed issue in the theory of high- 
temperature superconductivity is the question whether 
strongly correlated electron models in 2D share the most 
striking feature of their II? counterparts, spin charge sep- 
aration. Whereas in a Fermi liquid the low lying spin 
and density excitations can be modelled by particle-hole 
(i.e. composite) excitations of the fermionic quasiparti- 
cles, which represent the 'most elementary' excitations 
of the system, the situation is reversed in a spin-charge 
separated system. There, the soliton-like spin and den- 
sity excitations are themselves the elementary excita- 
tions, and the physical electrons can be considered as 
composite excitations. 

By numerical diagonalization one can obtain the exact 
excitation spectra of small clusters and should thus, at 
least in principle, be able to decide whether they are con- 
sistent with any given theoretical scenario. We therefore 
have performed a study of the dynamical spin and den- 
sity correlation functions of the t— J model, which reads 
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The Si are the electronic spin operators, c\ a — c| CT (1 — 
rii _ CT ) and the sum over < i,j > stands for a summa- 
tion over all pairs of nearest neighbors. All results to 
be presented below have been obtained for ID and 2D 
16 site clusters, with the parameter value J/t=0A. In 
the present study we restrict ourselves to intermediate 
and low electron densities, where the results are more 
easily understood than near half-filling [fil. To be more 
precise, in this doping regime the excitation spectrum of 
the 2D t—J model turns out to be completely consis- 
tent with the Fermi-liquid derived particle-hole picture, 
even deviations from the noninteracting case are small. 



Our results suggest that throughout the range of dop- 
ings studied, the t—J model represents a not even very 
strongly correlated Fermi liquid. No detail of our results 
would necessitate the introduction of the concept of spin 
charge separation. 

Using the standard Lanczos algorithm we computed the 
dynamical spin (SCF) and density (DCF) correlation 
functions: 
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C a (q,u) = -9f(* |Ot 



1 



w-(H-E )- i0+ 



Here j^o) (Eo) denotes the ground state wave function 
(ground state energy), and as the operator O we choose 
the Fourier transform of either the density operator n< f + 
n i,l ( a = or the spin operator (1 — p e )(m,-\ — fii,{) 
(a = s). We introduced an extra prefactor of 2(1 — p e ) 
for the z-spin operator because with this definition both 
correlation functions obey the same sum-rule, 

roc 

E/ dhjC a {q,u) = p e (l - p e ), 

which faciliates their comparison. In a system of nonin- 
teracting particles, the above two particle Green's func- 
tions can be expressed as the convolution of the single 
particle photoemission (PES) and inverse photemission 
(IPES) spectrum, i.e. 

C c (q,uj) = E I ^' A PEs{k,u3')A IP Es{q-k,u)-u'). 

k J 

Having computed the PES and IPES spectra by the 
Lanczos method we can replace the ^-peaks in these func- 
tions by Lorentzians and obtain the convolution numeri- 
cally. Thereby the brodening of the <5-peaks in the spec- 
tral function was taken 1/2 of the one for the correlation 
functions, because the numerical convolution effectively 
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doubles the width of the peaks. The convolution obeys 
the sum rule 




where n(k) = (c k a Ck,<j} denotes the ground state mo- 
mentum distribution and m(6) = (1 + 8)/2 with 5 the 
hole concentration. The last term is different from but 
turns out to be quite small, so that all spectral func- 
tions under consideration have a comparable integrated 
weight. 

Approximating the spin and density correlation function 
by a mere convolution is a drastic approximation, which 
neglects for example all vertex corrections; as we will 
show now, for the electron densities under consideration 
it is nevertheless a remarkably good approximation. Figs. 

compare the convolution estimate with the exact SCF 
and DCF for electron densities between 10/16 and 4/16. 
While the agreement for 6 holes is not yet very impres- 
sive (particularly for the SCF), the agreement becomes 
better and better with decreasing electron density. With 
the exception of the low energy peak structure at (tt, 0) 
and (tt/2, 7r/2) (which is determined by symmetry-related 
selection rules unaccounted for in the convolution, see be- 
low) the overall shape of the exact correlation functions is 
reproduced quite well, sometimes to the degree of a one- 
to-one correspondence between the dominant features. 
While SCF and DCF are not obviously identical, they 
both can be viewed as being derived by slight but differ- 
ent modifications of the convolution. 
We now turn to the low energy peaks at (n, 0) and 
(tt/2, tt/2), indicated by arrows in Figs. |||4|. Obviously 
these are not reproduced well by the convolution. We will 
show that this is the consequence of point group selection 
rules, which are not accounted for in the convolution. As- 
suming the existence of a free electron like Fermi surface, 
the 6 hole ground state would correspond to a closed-shell 
configuration, where the momenta (0,0), (±7r/2, 0) and 
(0,±p«/2) are occupied. The ground states with lower 
electron densities then would have 'holes' in the outer 
shell of momenta, which are therefore only partially oc- 
cupied. This is confirmed by the spectral function at 
(tt/2, 0), which shows only a PES peak immediately be- 
low the Fermi energy for 10 electrons, but both a PES 
peak immediately below and an IPES peak immediately 
above the Fermi energy for 8, 6 and 4 electrons. The low 
energy peaks in the convolution for momentum transfer 
(tt, 0) and (7r/2, tt/2) then originate precisely from tran- 
sitions within this shell of partially occupied momenta 
(±7r/2,0), (0, ±7r/2), as indicated in Fig. |a. We adopt 
the hypothesis that the lowest states of the system can 
be described by a 'cluster version' of Fermi liquid the- 
ory. To be more precise, we assume that e.g. for the 8 
electron case there is a one-to-one mapping of e.g. the 
lowest states of the cluster and 'quasiparticle states' of 
the type au^ak 1 ,i\FS) , where ak, a annihilates a quasipar- 
ticle, | FS) denotes the closed shell Fermi sea of 10 non- 



interacting electrons and k and k' are restricted to the 
shell of momenta (±7r/2,0), (0,±7r/2). We next assume 
that there is a residual interaction between the quasi- 
particles, which lifts the degeneracy between the various 
states of this type with given total momentum. To sim- 
plify the notation we label the momenta as indicated in 
Fig. ||b. Then, from the fact that the exact 8 electron 
ground state is a spin singlet and has d x 2_ y 2 symmetry, 
we conclude that it should be modelled by the state 

^( a i,T a 3,l - a U a 3,T - a2,ta4,i + a 2 ^a 4 . 1 )\FS) . (1) 

This would also be consistent with the pronounced d x 2_ y 2 
pairing correlations found for larger values of J/t at this 
electron density Next, the only states in our sub- 

space with momentum (tt/2, tt/2) are 

-^{a3^a 4 ,i±a 3d a 4A )\FS). (2) 

The antisymmetric combination is a spin singlet and even 
under reflection by the (1, 1) direction, so that its point 
group symmetry is incompatible with that of the ground 
state: there can be no peak in the DCF. The symmetric 
combination is a spin triplet and has an acceptable point 
group symmetry so that we expect a peak in the SCF. 
Both predictions are consistent with the exact spectra. 
The only possible state with momentum (tt, 0) and the 
required point group symmetry (even under reflections 
by both x and y axis) in the model space is 

-^{a 1A ai^ + a 3A a 3A )\FS). (3) 

This is a spin singlet so that the SCF cannot have a low 
energy peak. Since the point-group symmetry of this 
state is consistent with the selection rules, we expect a 
low-energy peak in the DCF. Again both predictions are 
consistent with the exact spectra, although the weight of 
the low energy peak in the DCF is extremely small. 
Next, literally the same explanation holds for the case of 
4 electrons, if we replace the annihilation operators by 
creation operators, choose the fully occupied T-point as 
the 'Fermi sea' and replace (3,4)^(1,2) in (||). This is 
because like the 8 electron ground state, the 4 electron 
ground state is a spin singlet with d x 2_ y 2 symmetry. We 
thus expect a low energy peak in the DCF at (tt, 0) and 
in the SCF at (tt/2, tt/2); inspection of the exact spectra 
shows, that this is indeed realized. 

We turn to the case of 6 electrons. Here we have to 
use quasiparticle states with 4 holes in the 10-electron 
closed shell. The exact ground state has unusual quan- 
tum numbers: its total spin is S — 2 and it belongs to the 
i?2 (or d xy ) representation of the C 4v point group. This 
suggests to model the S z — ~2 member of the ground 
state multiplet by the state a±,-[G2, -[as^a^lF S), which 
obviously has the correct point group symmetry. The 
S z = member of the ground state multiplet (which we 
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are considering) can now in principle be obtained by act- 
ing twice with the spin raising operator. It is easy to 
see, however, that all necessary conclusions can already 
be drawn from the final states: introducig the opera- 
tor ti j — l/\/2(ai, jo/j + a^a^i), which creates a spin 
triplet on the momenta i and j, the only state with total 
momentum (n, 0) and the required point group symmetry 
(odd under reflection by both x- and y-axis) reads 



-i=(ai, T ai,j. - az^a3,i)t2A\FS). 



(4) 



This is a spin triplet, so that we expect a peak in the 
SCF but not in the DCF. Similarly, the state with max- 
imal spin and momentum (7r/2,7r/2) that is even under 
reflection by the (1,1) direction reads 



(5) 



Again, this is a triplet, so that we expect a peak in the 
SCF but not in the DCF. All predictions are consistent 
with the numerical spectra. 

We thus have shown that the low energy peak structure 
can be understood completely by a simple 'cluster Fermi 
liquid theory', which relies on nothing more than ele- 
mentary selection rules. To further strengthen the evi- 
dence for this interpretation, we now compare the elec- 
tron momentum distribution n(k) = (ct a Ck,a) m the 
ground state (GS) to that of the final states (FS) asso- 
ciated with the low-energy peaks in the SCF and DCF, 
i.e. we consider 



Ara(fc) = n FS (k) - n GS (k). 



(6) 



Based on the above 'model wave functions', we can pre- 
dict the changes of the quasiparticle occupation numbers 
which accompany the respective spin or density excita- 
tion. These are listed in Table |. Up to a constant factor 
due to the wave function renormalization constant Z, the 
exact An(fc) then should be the same. 
To obtain the final state wave functions, we take ad- 
vantage of the fact that the energies of the low energy 
peaks give highly precise estimates for the eigenenergies 
E final of the SCF or DCF final states. Thus, by apply- 
ing powers of (H — E final)" 1 to some randomly chosen 
trial state (thereby employing the conjugate gradient al- 
gorithm) we can efficiently converge out the respective 
final state wave function []4| and obtain its momentum 
distribution. Then, Figs. 7j-|] show the An(k) for all 
low lying peaks in both SCF and DCF. By comparison 
with Table | it can be seen that the exact results are in- 
deed completely consistent with the particle-hole picture: 
the changes are always substantially larger for the Fermi 
momenta than for any other momentum, the reduction 
of the magnitude of An(k) as compared to the 'free elec- 
tron' estimate suggests values of the quasiparticle weight 
Z between 0.6 (for 8 electrons) and 0.8 (for 4 electrons). 
The losses of n(k) at the Fermi momenta always exceed 



the gains; this is probably due to an enhanced scattering 
in the final (=excited) states. We thus have rather un- 
ambiguous evidence for the validity of the Fermi liquid 
like particle-hole picture at low frequencies. For higher 
frequency the good agreement between the exact correla- 
tion functions and the convolution strongly suggests the 
same. 

We now turn to a comparison with the ID model, 
where spin charge separation rather than Fermi liquid 
behaviour is known to be realized at all doping levels || . 
Fig. ^| compares the SCF and DCF with the convolu- 
tion of the single particle spectral function for p e =0.5. 
The Fermi momentum is fci?=7r/4, and SCF and DCF 
show their lowest energy excitation at different momenta, 
2fci?=7r/2 and 4fci?=7r, respectively, which is characteris- 
tic of spin charge separation ||. Unlike in 2D neither 
correlation function is approximated well by the convo- 
lution. The latter consists almost entirely of structureless 
continua, whithout the sharp peaks of the exact correla- 
tion functions. 

For the low density regime our results are consistent with 
those of Jagla et al. || . These authors studied the group 
velocity of spin and charge excitations in small clusters 
and found different velocities for spin and charge in ID 
(indicating spin charge separation) but identical veloci- 
ties (indicating absence of spin charge separation) in 2D. 
On the other hand, evidence for spin charge separation 
at electron densities 0.75 and 0.2 was seen by Putikka 
et al. in the static SCF and DCF of the 2D t-J 
model obtained by high-temperature expansion. How- 
ever, a subsequent cluster study of the same quantities 
by Chen et al. || could not strengthen this evidence: the 
DCF showed no indication of a singularity at the 'char- 
acteristic wavevector' for spinless Fermions. Moreover 
the dynamical SCF and DCF differ drastically for low 
and high doping H so that it seems not very plausible 
that they both can be described by the same physics. 
Our results show that while for low electron density the 
peak structure of SCF and DCF is similar and consistent 
with the convolution, there arc differences in the pole 
strengths. Thus, while SCF and DCF apparently can be 
explained by particle-hole transitions between analogous 
states, the matrix elements are different. This would sup- 
port an explanation of the the static correlation functions 
e.g. by RPA calculations for the large-C/ Hubbard Model 

I- 

In summary, we have studied the spin and density corela- 
tion function for intermediate and low electron densities 
in the t—J model. For a spin-charge separated system, 
(i.e. in ID), these two correlation functions cannot ade- 
quately be described in a particle-hole excitation picture. 
Contrary to this, our results for 2D show that the sim- 
plest estimate based on the particle-hole excitation pic- 
ture, namely the convolution of the single particle spec- 
tral function, provides a remarkably good description of 
these correlation functions. The low energy peak struc- 
ture can be completely explained by a simple 'cluster 
Fermi liquid theory' which relies only on the particle- 
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hole picture and elementary selection rules. The respec- 
tive particle-hole transitions can be directly made visible 
in the momentum distribution. Our exact results thus 
suggest that for the less than quarter filled case the t—J 
model has a fairly conventional particle-hole type excita- 
tion spectrum, and represents a not even very strongly 
correlated Fermi liquid. There is no need to invoke spin- 
charge separation to understand the dynamics of the 2D 
t—J model for the electron densities under consideration. 
It is interesting to note that for electron densities near 
half-filling (i.e. 2 and 4 holes in the 4x4 cluster) the 
situation changes completely: the DCF consists almost 
entirely of incoherent, high-energy continua, the SCF has 
sharp low energy peaks and almost no continua . The 
dominant spin excitation at (w, ir) is a spin-wave like col- 
lective mode, i.e. a remnant of the undoped system; only 
spin excitations with different momentum transfer still 
correspond to particle- hole transitions p0| . In fact, quite 
a number of physical quantities show a significant change 
at hole concentration ~0.2-0.3: the doping dependence 
of the Drude weight fU]| , the temperature dependence 
of the susceptibility |12| |, the sign of the Hall constant 
p3| . For J=0 Chiappe et al. have shown that the 
overlap of the exact cluster ground state with Gutzwiller 
wave functions drops sharply from ~1 for low and inter- 
mediate electron densities to ~0 for high density. The 
clear distinction which can be seen in a wide variety of 
physical quantities suggests to assume two quite differ- 
ent phases for the 2D t—J model as a function of doping: 
for p e < 0.7—0.8 there seems to be a fairly conventional 
Fermi liquid, with a ground state that is 'adiabatically 
connected' to the nonintcracting one and correlations are 
of little importance. Near half-filling, the correlations 
dominate and the system seems to be continuous with 
the insulator: the spectral function shows rigid-band be- 
haviour pl l with increasing doping , the Fermi surface 
takes the form of hole pockets 0,[17| , the spin excitation 
spectrum is reminiscent of the undoped Heisenberg anti- 
ferromagnet |To| . 
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FIG. 1. Spin correlation function (right column, dotted 
line) and density correlation function (left column, dotted 
line) as compared to the convolution of the single particle 
spectral function (full line). All spectra are for the 4x4 
cluster t—J model with 10 electrons and J ft — 0.4. 

FIG. 2. Same as Fig. [i] for 8 electrons. 



FIG. 4. Same as Fig. |lj for 4 electrons. 

FIG. 5. Possible particle-hole transitions between the 
'Fermi momenta' in the 4x4 cluster (a), labelling of the 
Fermi momenta' (b). 



FIG. 6. Difference of n(k) for the final states belonging to 
the low energy peaks indicated by arrows in Fig. ^ and n(k) 
fort the ground state i.e. the quantity An(fc) defined in M). 
The figure shows all momenta in the Brillouin zone, the 'Fermi 
momenta' of Fig. |B| are indicated by dark symbols, the top 
right corner corresponds to (7r,7r). The number of electrons 
is 8. 



FIG. 7. Same as Fig. |6| for 6 electrons. 
FIG. 8. Same as Fig. |^ for 4 electrons. 

FIG. 9. a: Comparison of convolution (full line) and exact 
DCF (dotted line) for the ID 16-site chain with 8 electrons 
and J/t = 0.4. 

b: same as (a) but for the SCF. 



TABLE I. Change of the electron occupation numbers of 
the Fermi momenta as obtained from the model wave func- 
tions for the SCF and DCF final states. 
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FIG. 3. Same as Fig. hi for 6 electrons. 
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Fig. 4 
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